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Abstract 

In this publication I discuss the phase diagram of a frustrated spin-1/2 Heisenberg model sug- 
gested in A. A. Nersesyan and A. M. Tsvelik, Phys. Rev. B67, 024422 (2003). The phase diagram 
contains (w, 0) and (jr, ir) antiferromagnetic phases separated by the Valence Bond Crystal (VBC) 
state. I argue that the point of the phase diagram with deconfined spinons, predicted in the afore- 
mentioned work, is situated in the middle of VBC state, at the point where the dimerization order 
parameter changes sign. 

PACS numbers: PACS No: 71.10.Pm, 72.80.Sk 
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I. INTRODUCTION 



In our previous paper we described a model of spin-1/2 antiferromagnet (the so-called 
Confederate Flag or CF model) where the fine tuning of interactions gives rise to a state with 

n n n 

Pactional quantum spin number excitations (spinons) The subsequent studies (|2J,[3J, 
^] and especially 0) have provided a support to our results. In the original publication 
we did not discuss what happens if one deviates from this special point. Such discussion 
is a subject of the present publication. Here I discuss the phase diagram of CF model 
concentrating primarily on the vicinity of the deconfinement (D) point. This gives us a 
better understaning of the physics involved and also helps to put the CF model in the 
broader context of studies of frustrated magnetism. Though the corresponding literature is 
enormous (see, for example, 6( for review), the theoretical efforts are primarily concentrated 
on idealized models (such as models of dimers or gauge field theory models) whose relation to 
microscopic models with realistic interactions is not clear. The standard argument invokes 
universality: the belief is that low energy behavior of such systems will be independent on 
microscopical details following some universal patterns. It is always interesting to check 
general considerations against concrete models. In that sense, CF model, being simply a 
model of a Heisenberg magnet with short range interactions, presents an almost unique 
example. 

The literature knows two scenario for realization of the D-point. In one of them deconfined 
spinons exist on the boundary between antiferromagnetic and VBC state in the hypothetical 
situation when these states touch each other at a Quantum Critical Point (tJ. In the other 
scenario, based on the study of models of quantum dimers, spinons appear at the boundary 
between two VBC states (the Roksar-Kivelson critical point) 0, jj. As I shall argue in this 
paper, neither situation is realized in CF model. In that model the D-point separates two 
VBC states, as in the second scenario, but the analogy does not go much further. First, 
according to [sj],[i3], the spectrum at the D-point in the dimer models consists of spinless 
particles ("photons") with a quadratic u ~ k 2 spectrum. This is absolutely incompatible 
with CF model which is approximately (1+1) Lorentz invariant when the interchain coupling 
is weak. This symmetry dictates that the dispersion of the gapless excitations along the 
chains direction must be linear. Second, one of the VBC phases in the dimer models contains 
a 'devil's staircase' of commensurate and incommensurate phases which does not agree with 
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the situation in CF model where both VBC phases are simple. So it appears that CF model 
is quite distinct and does not fit into the known categories. 



II. THE MODEL 

The model suggested in our original paper [lj is a spin-1/2 Heisenberg magnet with 
spatially anisotropic exchange interactions. The exchange in one direction is much stronger 
than in the others and therefore this model can be viewed as a model of weakly coupled 
chains. Recently Batista and Trugman p] have found an isotropic version of the CF model 
and shown that it has the same ground state degeneracy and possesses spin-1/2 excitations. 
Therefore the requirement of the space anisotropy is only a matter of convenience. It allows 
us to use the continuum limit in one direction and to apply the field theory methods. The 
existence of the deconfined point also does not depend on the number of transverse directions, 
therefore for the sake of simplicity I will discuss the two-dimensional version of the model. 
In that case the interaction pattern reminds the Confederate Flag (see Fig. 1). The CF 
model Hamiltonian is given by 

j,n ' &j+l,n + ^ ] [J r Sj tTl + Jd (Sj_|_i 5ri + Sj_i )n )] • Sj jn _|_^ / + Wquarten (1) 
j,n { fJ,=±l ) 

where Sj )n are spin-1/2 operators, and J\\ » J r , Jd > . The term V quarter contains a four- 
(and possibly higher) spin exchange interaction with a small coupling constant A ~ J? / J\\- 
This term was absent in the original publication, but, as was recently demonstrated by 



Q. 



one needs to introduce it to fine-tune the model to the state with 



Balents and Starykh 
deconfined spinons. 

Assuming that the interchain couplings (J n Jd, A) are much smaller than the exchange 
along the chains (Jn) it is legitimate to adopt a continuum description of individual chains. In 
this description, the local spin densities are represented as sums of the smooth and staggered 
parts: 

Sj,n/ao -> S n (x) = Mn(x) + (-iyN n (x), x = ja Q , (2) 

ao being the lattice spacing in the chain direction. 

The low-energy dynamics of the spin-1/2 Heisenberg antiferromagnet 

H id = J \\ H( S i S i+i) ( 3 ) 



FIG. 1: Exchange interactions pattern for CF model, the red lines correspond to J r and the green 
ones to Jd- 



FIG. 2: A possible realization of the Confederate Flag exchange pattern via a superexchange. The 
superexchange between spins of magnetic (black) ions occurs through orbitals of non-magnetic 
(green) ones with a big ionic radius. The pathes 1 and 1' generate the vertical interchain exchange 
( J r ) and the pathes of type 2 generate the diagonal one such that J g = J r /2. 

is described by the SUi(2) Wess-Zumino-Novikov-Witten model. The latter Hamiltonian 
can be written in terms of the so-called chiral vector current operators, J and J, satisfying 
the level k = 1 Kac-Moody algebra (this approach has been described in a vast number of 
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publications; (for a review see 



H or |ll | 
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ItXV f r — — 

H 1D ^—jdx [: (J- J) : + : (J- J) : 

with v = 7ra J||/2. 

It is remarkable that the smooth part of magnetization, 

M = J + J, 



+ ■ 



(4) 



and the spin current, 



v(3- 3) 



(5) 



(6) 



are locally expressed in terms of the chiral currents. 

In the CF model, the exchange is frustrated in the direction perpendicular to the chains. 
In order to eliminate the coupling of the leading relevant operators (the staggered mag- 
netizations N n (x) • N n+ i(x) and the dimerizations e n (x)e n+fl (x)) one has to fine tune the 
couplings J r — 2J d and A. In the leading order in J r , J d the fine tuning is achieved when 
J r - 2 J d = 0, A = 0. 

To discuss the phase diagram of CF model one needs to deviate from the fine-tuned point 
with deconfined spinons. For weak interchain interactions one still can employ the continuum 
Hamiltonian obtained using the continuum description of individual chains following the 
asymptotic representation (J2J) of the spin operators: 



H = H l + H 2 , 

2N 

n=\ 
1 

2 



z n=±i 



H- 



H\D,n 



(7) 
(8) 

(9) 



where Hm is given by Eq.(0]). For the lattice Hamiltonian presented on Fig. 1 we have 
7 = J r + 2Jd and SJ = J r — 2J d . One very important property of model (jBJ) is its (1+1)- 
dimensional Lorentz invariance which survives at 5J ^ 0. This dictates the form of the 
excitation spectrum for all particles: 



E 2 = v 2 k 2 + A 2 (k A 



(10) 



where A(kj_) is a periodic function of the wave vector component perpendicular to the chains. 
I emphasise that this result is valid for an arbitrary number of chains. As I have mentioned 



in Introduction, this property precludes the existence of k 2 gapless modes characteristic for 
dimer model critical points. 

Model (|8I9|) can be viewed as the critical system (the bunch of non-interacting spin- 1/2 
chains) perturbed by relevant interactions. Each of these interactions generate their own 
energy scale according to their scaling dimension. The current-current interaction, being 
only marginally relevant, generates the scale 

M ~ J|| exp[-7r 2 J\\/2{J r + 2J d )\ (11) 

and the interaction of the staggered spin components generates a scale ~ 8 J. When \8J\ » 
M the system orders antiferromagnetically (though in two dimensions only at T = 0). The 
corresponding Neel wave vectors are (ir, 0) for 5 J < and (tt, tt) for 5 J > 0. These things are 
absolutely obvious; it is less obvious however what happens in the opposite limit M » \5J\. 
For weak interchain interactions the best I can do in this case is to study the four chain 
model, where non-perturbative calculations can be carried out explicitely. I believe that the 
four chains are representative enough to give an insight into what happens for an infinite 
system. 



III. THE FOUR CHAIN MODEL AS BOTH SOLVABLE AND REPRESENTA- 
TIVE CASE. 

So let us consider the case of four chains with periodic boundary conditions in the trans- 
verse direction. Let me briefly recall the results for 5 J = obtained in As was 
noticed that at 5 J = the model acquires an additional symmetry: the sectors with differ- 
ent parity decouple. This follows from the fact that the relevant interactions couple only 
currents of different chirality belonging to different chains. Thus the relevant current-current 
interaction is 

(Ji + J 3 )(J 2 + J 4 ) + (J 2 + J 4 )(Ji + J 3 ) (12) 

and the Hamiltonian (JHJ) decouples into two parts: 

H = H + + H (13) 

The + parity sector contains Ji j3 and J 2j 4 currents and the — parity sector contains Ji j3 
and J2,4 currents. Each of the models represented by the Hamiltonians H± is integrable, 
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which can be demonstrated by rewriting them in terms of familiar integrable models. Let 



us recall how it was done in 



moreover, according to the sum of two SUi (2) WZNW models (the central charge 2) can 



A sum of two k — 1 SU(2) currents is the k = 2 current; 
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be represented as the SU2(2) WZNW model with central charge 3/2 and plus one mass 
Majorana fermion (a critical Ising model) with central charge 1/2. Using the results of 
we rewrite the entire Hamiltonian density (jlj) as follows (here only the (+)-parity part is 
written) : 

+ H massive 
if if 

Hmassless = —^XO^xXO + ^XQ^xXQ (14) 
77*71 _ _ _ 17; fy 

H massive = — (: I • I : + : I • I :) + 7 I • I = -(- X a d x X a + Td.T) ~ \{x a X a ? (15) 
where a = 1, 2, 3 and 

I = Ji + J 3 , I = J 2 + J 4 

Ji,3 = ~{±XoX + gftx Xl}» ^2,4 = ^ |±XoX + \[X x X]} (16) 

The fields x, x stand for real (Majorana) fermions. 

Eq.((T3J) describes a critical Ising model; the corresponding excitations are gapless and 
non-magnetic; they appear in the sectors with both parities. 

Let me say several words about the 0(3) Gross-Neveau model fjl5|) . Though Majorana 
fermion description presents some advantages, the staggered magnetization components as 
well as their product Q are nonlocal with respect to these fermions. It turns out that 
the latter interaction can be expressed in terms of order and disorder parameter operators 
of the eight Ising models corresponding to each Majorana fermion species. Since we are 
interested in the case when \8J\ < M, we have to recast the perturbation in terms of the 
nonchiral fields of models H±. This can be done using Abelian bosonization representation 
for individual chains and the correspondence between C = 1 theory and two critical Ising 
models (see [llj|, The net result for the density of the perturbation Hamiltonain 

obtained after some algebra is 



H 2 = 5J(N 1 + N 3 )(N 2 + N 4 ) = 

5J( y fi^fi )[3(a 1 a 2 cr 3 ) + ( y aia2a 3 y + {^2^Y " (^1^2/^3)" + •••] (17) 
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where the dots stand for the terms which do not have finite averages at 5 J. ut,^t ( a = 
0,1,2,3 are order and disorder parameter operators of the Ising models associated with 
the Majorana fermions Xai Xa from + and — sectors. As we shall demonstrate later, this 
interaction leads to confinement of spin-1/2 particles. 
The dimerization order parameter is 

(Ni-N 3 )(N 2 -N 4 ) = 

( cr o" (T o")[( cr l (T 20-3) + (o"lO'20'3)~ + 3(/ii/i 2 /i 3 ) + (/il/U2/i3)~ + •••] + 

04 A*0 )[-3( (T l "20"3) + (0-l0"20-3) _ + (^2^) + (j^lfMlUs) ~ + •••] (18) 

As I have said, the terms in the square brackets in Eq.(|17|) have nonzero vacuum averages 
even at 5 J = 0. The unperturbed model has four ground states: with {(<Ji<J2<J3) ± ) ^ 
(the a ± vacua) and ((/ii/x 2 A i 3) ± ) ^ (the ^ vacua). Replacing in (fTTj) the corre- 
sponding products by their vacuum expectation values ~ M 3//8 and identifying (ixqUq) = 
cos( v /7T$ ), (ao (Tq) = sm(^^ ), we get 

H 2 ~ (5J)M 3/4 costv^F^o + tan" 1 3] (( / u ± ) ^ 0) 

~ (5J)M 3/4 cos[v^F$o - tan" 1 1/3] ((a ± ) ^ 0) (19) 

and if for one parity (a) ^ and (//) ^ for the other. Thus the ground state degeneracy 
is now reduced to two; the order (disorder) parameters in both sectors now condense simul- 
taneously. The o — > ji degeneracy is not lifted, since two forms of the potential ((EH) are 
equivalent under the uniform shift y^^o ~~ * V^^o + tt/2 which does not affect the gradient 
term (d^o) 2 . Perturbation Q couples the critical Ising model sectors with different parity; 
being projected on the state with ^ the resulting Hamiltonian becomes the 

sine-Gordon one: 

H SG = \xt^d xX t + ^Xor 3 d xX o + AM^(5J)(a+a + 3/i Vo ) = 

Uu 2 + (d x %) 2 ] + A(5J)M 3/4 cos^^o + tan" 1 3) (20) 

where A ~ 1 is a numerical coefficient, r 3 is the Pauli matrix and n is the momentum 
density operator. The cosine term gives rise to the mass gap 

A ~ (<5J) 4 / 7 M 3 / 7 (21) 
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The average dimerization is 



e = ((N 1 -N 3 )(N 2 -N 4 )> = 

~ M 3 / 4 (cos[V^$o - tan" 1 3]) ((//) ^ 0) 

~ -M 3/4 (cos[0r$o + tan" 1 1/3]) ((a) ^ 0) (22) 

This average does not depend on the choice of vacuum (that is whether o or fi are in the 
condensate) : 

e ~ (6J) 1/7 M 6 / 7 (23) 

Thus the dimerization changes its sign when S J goes through zero following a power law 
dependence. The exponent is quite small and it is possible that in the limit of infinite 
number of chains the transition is the 1st order, as in the isotropic CF model Q]. I believe 
that this is indicative of the physics behind the deconfinement: it occurs on the boundary 
between two Valence Bond crystalline orders. 



IV. CONFINEMENT OF SPINONS 

The problem of kinks confinement in the potential (fT7|) is somewhat peculiar differing 



3, Q, Q 



. In our 



from the standard confinement problems studied by various authors 
case the string tension is provided by the cos[-y/7r$oL sin[-y/7r<I>o] terms and becomes energy 
dependent. Though it is not directly related to the main topic of the paper, it is interesting 
enough on its own right. So I will spent some time discussing this problem. 

Let us consider a kink interpolating between the a- and fi- vacua in the + sector ([/icx] + 
kink). It carries an isotopic (spin) index a = ±1/2. According to (J17)) . the creation of such 
a kink will lead to the rise of the total energy proportional to the system size. To prevent 
this, one has to create a similar kink in the — sector. If we have two kinks fia - one centered 
at x — x\ and the other one at x = x 2 , the effective potential for the = y^^o + tan _1 3 
field is 

A(5J)M 3/4 {e(x x - x) cos[v^F$] + 9(x - x 2 ) sin[v^F$]} (24) 

The kinks are heavy particles whose mass M far exceeds the masses of excitations of the 
$ field ~ A. This justifies the approximation which takes the kink configuarations as step 




FIG. 3: The dimerization patterns for two VBC orders. The dots are spins; the strong bonds are 
depicted in color. 

functions. At x < X\ field $ is locked at y/rr, at x > x 2 it is locked at 3y/n/2. The solution 
in the middle is 

i~~ 7T X — X\ 

\/7r$ 



2 x 2 — X\ 

such that the total energy difference between the vacuum without kinks and the vacuum in 
the presence of two kinks is 

V = B 2 A 2 \x 12 \ + -^- (25) 
8|xi 2 | 

where B ~ 1 (though its numerical value can be calculated, it is not of an interest here) and 

B 2 A 2 = A(5J)M 3 ^(cos[^]) (26) 

This energy difference provides an effective confining double-well potential in which the 
[/xcr] + — [[ia}~ bound states are formed. Recall that the potential does not depend on the 
spin configuration (this will be no longer the case when A ~ M). Therefore the bound 
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FIG. 4: The schematic phase diagram for the CF model. The bold green lines denote the 1st 
order phase transitions between the antiferromagnetic and the dimerized states, x = 5J/M. The 
dimerization changes sign at x = 0. 




X 



FIG. 5: The schematic form of the confining potential (|25|) . The energy levels corresponding to 
odd and even parity states are colored in blue and red respectively. 

states form 4-fold degenerate isotopic multiplets consisting of SU(2) triplets and singlets. 
The mass spectrum of kink bound states in such double-well potential is somewhat different 
from the purely linear confining potential studied in jlfi| . \v\ . jisl] : 

M n -2M & 
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5A 1^/2+ L/s/tt 



(n + l/2)[l±<J(n)] + ...L ri«(M/A) 2 / 5 (27) 



.1/ 



where 



~ exp -(tt/2) 17 / 8 (M/SA) 3 / 4 



1 



n+1/2 



and 



M n - 2M » BA(n 2 BA/M) 



1/3 



n » (M/A) 2/5 



(28) 



Not all these bound states are stable; the particles with masses greater than 2M + A can 
decay into particles with smaller masses emitting excitations of the $o field. Therefore 



particles with n > yM/A are unstable. Since the power 1/2 is quite close to 2/5, the mass 
sequence of Eq. (J28|) is never reached. 

As far as kink-antikink states are concerned, for them the string potential is not repulsive, 
but attractive and they are expelled from the spectrum. 

V. CONCLUSION 

The study of the four chain case gives reasons to believe that the infinite system has a 
phase diagram presented on Fig. 4. The deconfinement point appears in the middle of the 
VBC phase, as in the dimer models considered in j^,!^. The excitation spectra are different, 
however. The spinless modes of the dimer models have a quadratic spectrum uj ~ k 2 ; such 
spectrum cannot emerge in CF model due to the (l+l)-dimensional Lorentz invariance. 
Thus I conclude that the D-point of CF model belongs to a universality class different from 
the universality class of the Roksar-Kivelson critical point. 

In conclusion to this paper I would like to point out a rather curious parallel between the 
problem of frustrated magnetism and another long standing problem of condensed matter 
physics, namely the problem of heavy fermion state formation in rare earth compounds. In 
these compounds magnetic moments of highly localized electrons belonging to rare earth 
ions / shells interact with delocalized electrons from the broad conduction band. The prob- 
lem is discussed in the literature in terms of competition between the Kondo screening and 
the induced interspin interaction (the so-called RKKY interaction). On the formal level 
this competiton looks exactly like the competition between the less relevant current-current 
interaction in model (jHJ) and the more relevant interaction of staggered magnetisations Q . 
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Indeed, the energy scale generated by the Kondo screening (Kondo temperature) is exponen- 
tially small in the coupling constant as in (jllj) and the energy scale generated by the RKKY 
interaction is proportional to the square of the spin-fermion coupling. Though one would ex- 
pect that exp(— 1/g) is always much smaller than g 2 , there is a vast class of materials where 
the Kondo screening is manifest at temperatures much larger than the temperature of mag- 
netic ordering. One possible explanation is that the RKKY interaction with its oscillatory 
behavior, is highly frustrated which leads to cancellations similar to the one considered in 
this paper. 

I am grateful to O. Starykh, I. Zaliznyak, A. Nersesyan and F. H. L. Essler for discussions 
and interest to the work. I acknowledge the support from US DOE under contract number 
DE-AC02 -98 CH 10886. 
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